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The results of calculations of the spatially-resolved density of states (DoS) in an S(F/N) bilayer 
are presented (S is a superconductor, F is a metallic ferromagnet, N is a normal metal) within 
quasiclassical theory in the dirty limit. Analytical solutions are obtained in the case of thin F, 
N layers which demonstrate the peculiar features of DoS in this system. The dependencies of the 
minigap and the DoS peak positions on the exchange energy and parameters of the layers are studied 
numerically. 
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In the past few years there was a noticeable interest to 
the Josephson junctions with ferromagnetic barriers due 
to possibility to realize the tt— junctions having the phase 
difference tt in the ground state. The tt— states in SFS 
junctions were first predicted by2i^ and realized exper- 
imentally by Ryazanov et almt^ in Nb/CuNi/Nb struc- 
tures and later by other groups^&SiiSiii using different 
ferromagnetic barriers. These experiments stimulated 
further theoretical activity (sec"'^ for the review). In par- 
ticular, Josephson structures composed from arrays of 0— 
and TT— Josephson junctions should exhibit extraordinary 
characteristics Such arrays were recently realized in 
zig-zag HTS/LTS structuresii. 

The purpose of the present paper is to study spatially 
resolved electronic density of states (DoS) in the struc- 
ture of S(FN)type, consisting of a bulk superconductor 
with ferromagnetic and normal layers on the top of it, 
which is a generic system for 0- and tt— junctions con- 
nected in parallel. 

DoS in SF bilayers (a ferromagnet coupled to a 
superconductor) was studied quite extensively before. 
Two new features were predicted compared to SN sys- 
tems: spin splitting and spatial oscillations of DoS in a 
ferromagne1*iSiiLiSii2iSa*2iiS^. The effect of spatial oscil- 
lations was quite extensively discussed in the theoreti- 
cal literature in different modelsiSiiLiSiiS and observed 
experimentally^^. This effect is closely related to O-tt 
transitions. The effect of splitting is relevant for thin fer- 
romagnetic layers and was studied theoretically in2Si2i. 
In the present work we discuss an interplay between 
the oscillations and splitting in a more complex S(F/N) 
structure. 

The geometry of the structure is shown on Fig.l. We 
assume that the dirty limit conditions are fulfilled in all 
metals, F is a weak monodomain ferromagnet with the 
exchange energy H much smaller than the Fermi energy 
and the interfaces are not magnetically active. In this 
case spin dependent corrections to the resistivities can 
be neglected and the S(F/N) structure is described by 
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FIG. 1: The geometry of the structure. 

the following spin independent suppression parameters 
IBF ^ RbfAbf / Pf^f, If = PsS,s/ Pf£.f, (1) 



iBN = RbnAbn I PnS,n , In = PsCs/pnCn, (2) 



"IB = RbAb/ Pn£,Ni 7 = Pf^f/pn£,n- (3) 

Here Rbf, Rbn, Rb are the specific resistivities of the 
SF, SN and NF interfaces respectively; ps,f,n , Ds,f,n 
and S,s,F,N are the resistivities, the diffusion constants 
and the coherence lengths of the S, F and N layer and 
the coherence lengths, where £,s.f.n ~ \/ Ds.f,n/2ttTc 
and Tc is the critical temperature of the superconductor. 

Under the above assumptions the problem can be 
solved in the framework of the Usadel equations^^. To 
simplify it further we assume that S is a bulk supercon- 
ductor and 7Ar ^ 7BAr, 7f ^ 7b_f so that the rigid 
boundary conditions 



2 



are valid for superconductor. Here A is the magnitude 
of the order parameter in S electrode, Fs and Gg are the 
Green's functions, uj — ■KT{2n + 1) are the Matsubara 
frequencies. 

Let us choose the x^y axes as shown in Fig.l and use 
the 9 parametrization G ~ cos 9, F — sm9, then the 
Usadel equations have the form 



F,N' 



N 



dy- 



■9f. 



N 



sin6'F,Ar = 0, (4) 



where uj — uj -\- iH in F and cu ~ lo in N. 

The boundary conditions at the SF {y = 0, — cxd < a; < 
0), SN (y = 0, < X < oo) and FN (a; 0, < y < dp, 
(In) interfaces have the form^ 



the decay lengths and (f are 



^N,F 



^N,F\ 



' ttTcJ cos 0Ar,_F(±Oo) 
{uj J + ttTc COS 9s) 



(11) 



and we have taken for simplicity equal barrier parameters 
for F and N 



dN 

iBN-r- 



dp 

IBF-^ 



T (12) 
The general solution of Eq. @ has the form 

9nAx) = 9nA^^)+ (13) 



1b(f,n)^f,n-q;^9f.,n = - sin(6'5 - 9f,n), y = 0, (5) 



1b£,f-q^9f = sin(6'Ar - 9f), x 0, <y < dp, ^at. 



(6) 



d d 
^N^0N = 1^ft-9f, a; = 0, < y < dp^d^, 
ox ox 



where sin^s = A/Vw^ + and cos 9s = uj/^/uj'^ + A^. 
At the free interfaces the boundary conditions are 



dy 



'f,n = 0, y = dF,N 



(7) 



—0F,N--O, 0<y<dF, x^T^ (8) 

We will consider the limit of thin F and N layers 
dF.N ^ ^F,N- In this case one can neglect both the 
derivative on x and nongradient items in Usadel equa- 
tions I0J and substitute the resulting solutions 



dp 



9F,N{x,y) = 9f,n{x) - KF,N{x)'^^^y + Kf.n{x)- ^ 

'iF,N '^^F,N 



y 



Kf.n{x) = 



F,N 



into boundary conditions ((Sj. Then the problem is re- 
duced to the one-dimensional equations for lateral varia- 
tions of F in the x-direction: 
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(n.f-q^On,f - sin(6'7v,F - 9n^f{±'X))) ^ 0, (9) 



where 



9n,f{'^oo) — arctan 



ttTc sin 9 s 
Loj + ttTc cos 9s ' 



(10) 



-4 arctan 



(tan — ■ ' ' '-) exp{T7 } 



N.F 



The integration constants 9n (0) and 9f (0) in ( I13|l have 
to be determined from the boundary conditions lO at 
a; = and can be found analytically in the limit of large 
transparency of the FN interface when 9{x) is continuous 
at at X — 



9n{0) = 9f{0) = 9{0). 

From it follows that the effective decay length 
in normal metal, (^n, is a real quantity and equals to 
Cn = S,N\/j for small lj and tends to Cn = £,n\/t^Tc/uj 
with UJ increase. The effective decay length ^f in fer- 
romagnet and at low w <C A, H/j is given by (^f = 

^F\Jl/ — I'^iH/irTcY i. e. it becomes complex for 
sufficiently strong exchange field H > ttTc/j. 

Below we consider several limiting cases. 

Identical F and N metals. Assume for simplicity 
that the F and N materials differ by the existence the 
exchange field in F (7 = 1, = ^tv = C)i then from © 
for 9{0) we have 



9{0) = 2 arctan- 



sm ■ 



_««i^+gcos^^(^' ^~Cf 

(14) 

Identical F metals with antiparallel direction of 
magnetization. The results can be easily generalized to 
the case of an S(F/F) structure with two identical ferro- 
magnetic films having opposite magnetization directions 
(antifcrromagnetic configuration) 



9{0) = 2 arctan 



g sm 



g sm — - 



g* cos 



g cos ■ 



(15) 



Using the solutions obtained above one can calculate 
the spatially resolved DoS in S(F/N) and S(F/F) struc- 
tures. 
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FIG. 2: The total DoS in an SF bilayer for various values of 
7/1 as indicated in the figure. 



FIG. 3: Spin resolved DoS: comparison of FS, FN and S(FN) 
for h"f=0.5. 



DoS in S(F/N) and S(r/F) proximity systems. 

The DoS for each spin direction is given by 



Nis) 



No 



ReG{uj -ie + 6), 



(16) 



where Nq is the total DoS for both spins at the Fermi 
surface in the normal state and G{e — iS) = cos0(e — iS) 
is the retarded Green's function. The total DoS is found 
by summing over both spin projections i.e. ]V*°'°' = 
N{H) + N{-H). 

DoS in N and F metals far from the F/N in- 
terface. In a normal metal far from the F/N interface 
(a; = 00) the total DoS is given by 



NN{e) = NoRe 
£n - 



£(1 + 7VA2 



+ J^ 



P^~A^sgn{e) (17) 



It is well known (see Refsi ^^i^? ) that DoS in a F/N bilayer 
has a minigap at Eg < A, which depends on the value of 
7, and A^Ar(e) has the peaks at e = Eg and e = A. The 
minigap Eg characterizes the strength of superconducting 
correlations induced into N metal due to the proximity 
effect. 

In SF bilayers, modifications of DoS due to spin split- 
ting of energy levels were investigated in Refsi^SiSl. The 

DoS per spin projection in the F layer has the form 



^^F'\,Fi(£) - — He . 



(18) 



?2 



A2 sgn(e =F H) 



EF]',Fl = E + j{eT H)\/A^ ~E^, 



which demonstrates the energy renormalization due to 
the exchange field. In particular, it follows from l|18|) 



that now there are two minigaps in the spectrum £g| and 
Egl and egt <^g< ^gl- 

The total DoS in a F/N bilayer A^tot(e) ^Nfi{e) + 
Np^{e) is shown in Fig. 2. It is clearly seen that at 
h = H/ttTc < 1/7 there are three peaks in DoS lo- 
cated at Egi, Eg^ and A respectively. At /i = I/7 the low 
energy singularity is shifted to the Fermi level and for 
/i > 1/7 the first peak disappears resulting in only two 
singularities in the DoS at e = Eg^and e = A. Note that 
the total DoS at low energies depends nonmonotonously 
on H even in a thin F-layer, even though spatial os- 
cillations are absent across the layer. Eq. ((TH|) yields 
N^,,^,ie = 0) = iNo/2)Re^hsgnih)/y/WW^). For 
7ft, < 1 the total DoS N{0) = due to the minigap in 
F while for 7/1 > 1 the total low-energy DoS increases 
sharply, exceeds unity and saturates at N{Q) — Nq for 1 
jh > 1. 

DoS at the F/N interface. At x = Oand for iden- 
tical transport parameters on the F and N metals from 
((TE|l we obtain 



where 

£fn 



— Re — 
2 n 



-IEn — iSF'l^Fi + "^sfn 



FhFi 



N 



2v/A5 



^FN 



^Nflp^^Fl — EnEf^,FI 



A2 



(19) 



(20) 



and flNjEN and flFi,Fi,£F^,Fi are defined by IjlTfl and 
(|18|l respectively. 

It follows from Ea. (|19|l that, similar to the case of the 
SF bilayer considered above, the minigap exists if 7/1 < 1. 
With increasing exchange field the total DoS at e = be- 
comes nonzero if 7/1 > 1 at is given by simple expression 



A^(0) = NoVl^h^ - 1/7/1. 



(21) 



4 



SN 

SFup 

S(FN)up 

3- -o-S(FN) total 




J -| 1 1 1 1 1 1 T 1 1 1 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



FIG. 4; Spin resolved DoS: comparison of FS, FN and S(FN) 
for h'^—2. 



imaginary number and both numerator and denominator 
in Ea. l|19|) are complex thus leading to nonzero DoS in 
this energy range. Similar to the spin-down case, there 
is no peak in Npij(s) at the gap energy e = Eg, while the 
peak occurs at e = Sg^ (see Fig. 3). With further increase 




DoS at the F/F interface. At a; = we have 



— 2ieN + "^Sff fig. 5: The total DoS in S(FN) for various values of 7/1 as 



N{e) /No = Re — (22) indicated in the figure. 



eFF = sJ ^ 

It can be shown that DoS at the F/F interface given by 
Eq.(|221l coincides exactly with the total DoS for the F/N 
interface, Np^ie) + Npi{e), where Np^{e) and Npi{e) 
are given by Ea. (|19() . In particular, the minigap exists if 
7/1 < 1 and at 7ft, > 1 the DoS at F/F is determined by 
Eq.lini. 

The results of calculations from Ea. (|19|) at low tem- 
peratures T <^ Tc are shown in Figs. 3, 4 for = 0.5 
and ft,7 = 2, respectively, together with the DoS for SF 
(x —00) and SN {x — *■ 00) bilayers. 

There are four characteristic energies in the system: 
Egi, Eg, and A. Here £gi is the minigap for the spin- 
down subband SF bilayer at x — > — cx). It follows from 
Eq. (fT^ that Np^ (e) = at e < Sgi and becomes nonzero 
at e > £gi, i.e. e^j^ is the minigap for the spin-down 
subband in S(FN) at a; = 0. However, contrary to SF 
case Np^^s) has no peak e — Sgi but grows continuously 
from zero value. 

For the spin- up subband, the minigap in Np^ (e) is not 
equal to the gap Sgi in the spin-up subband in SF bilayer 
at X — > —00. Instead, Np^{e) the gap value is determined 
by Eg, the minigap in SN bilayer ai x 00. The formal 
reason is that in the interval e > Eg, flN becomes an 



of energy there is a sharp peak in DoS at e = A followed 
by saturation at No/2 for e 3> A. 

For /17 > 1 the minigap at Npi{e) vanishes and the 
structure of DoS becomes different, as illustrated in Fig. 4 
for the case — 2. The main qualitative difference from 
the previous case is that the spin-down and total DoS are 
gapless for /17 < 1. 

The total DoS at the F/N interface at a; = (which 
coincides with the total DoS in the F/F case), is shown in 
Fig. 5 for various values of ft,7.0ne can see that the gap is 
closed at ft,7 = 1 and the broad zero-energy DoS peak is 
formed with further increase of h until low-energy states 
become continuously filled at /17 ^ 1. 

In conclusion, we have studied theoretically the spa- 
tially resolved DoS in the S(FN) structures and in 
S(FF) structures with antiparallel magnetization direc- 
tions. Analytical solutions were obtained in the case of 
thin F, N layers which demonstrate the peculiar features 
of DoS in this system. We have illustrated the results nu- 
merically and have studied the dependencies of the mini- 
gap and the DoS peak positions on the exchange energy 
and parameters of the layers. 

This work has been supported in part by Russian Min- 
istry of Education and Science and RFBR Grant N 04 
0217397-a. 
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